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^— , Abstract 

Over the previous decade, many approaches for the modelling of radio emission from cosmic ray air showers have been developed. 
However, there remained significant deviations between the models, reaching from important qualitative differences such as unipo- 
Q-i/ar versus bipolar pulses to large variations in the predicted absolute amplitudes of up to factors of 20. Only recently, it has been 
^ realized that in the many models predicting unipolar pulses, a radio emission contribution due to the time-variation of the number 
of charged particles or, equivalently, the acceleration of the particles at the beginning and the end of their trajectories, had not been 
taken into account. We discuss here the nature of the underlying problem and demonstrate that by including the missing contribu- 

, tion in REAS3, the discrepancies are reconciled. Furthermore, we show a direct comparison of REAS3 and MGMR simulations for 

PjQ a set of prototype showers. The results of these two completely independent and very different modelling approaches show a good 
HpH level of agreement except for regions of parameter space where differences in the underlying air shower model become important. 
• This is the first time that two radio emission models show such close concordance, illustrating that the modelling of radio emission 
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from extensive air showers has indeed made a true breakthrough. 
Keywords: radio emission, extensive air showers, modelling and simulation 
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Due to renewed interest in radio detection of cosmic ray air 
showers, the need for reliable and detailed models of the radio 
emission process has increased continuously over recent years. 
Although the emission process had already been investigated 
extensively in the early days of air shower radio detection, it be- 
came clear that renewed modelling efforts were necessary. Con- 
sequently, modelling of the radio emission of EAS was restarted 
by a number of different groups as early as 2001. (A review on 
the evolution of the modern radio emission approaches can be 
found in fl].) 

When more and more models appeared on the market, it be- 
came obvious that there were contradictions between them. In 
particular, two classes of models could be identified. On the 
one hand there were models predicting unipolar radio pulses, 
the frequency spectra of which levelled off at a constant value 
towards very low frequencies. Most strikingly, this was the case 
for most of the models working in the time-domain. On the 
other hand, there were models predicting bipolar pulses, the 
frequency spectra of which fall to zero for very low frequen- 
cies. All frequency-domain approaches had this characteristic 
feature. 

These issues were first discussed in detail during a workshop 
dedicated to radio emission theory at Forschungszentrum Karl- 
sruhe H. In an article by Gousset, Lamblin & Valcares 
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a direct comparison was made between a description predict- 
ing bipolar pulses and a formulation producing unipolar pulses 
(see Fig. [T}. The resolution of the contradiction, however, was 
not yet found. 

From the fact that the source of the radio emission exists only 
over a finite time in a finite region of space, however, one can ar- 
gue that the zero-frequency component of the emission (which 
corresponds to an infinite time-scale), can contain no power 
fl. Therefore, the models with unipolar pulses, which have 
frequency spectra leveling off at a constant value at frequency 
zero, had to be suffering from some problem. The resolution of 
this problem has been found only recently. 

In this article, we shortly review the models producing unipo- 
lar and bipolar pulses before explaining what causes the dis- 
crepancy and how the affected time-domain models can be cor- 
rected. Finally, we compare two very different and completely 
independent models with correct implementations, namely the 
REAS3 JlH] and MGMR US] models. We demonstrate that 
for the first time, two models show a good agreement in the pre- 
dicted radio emission features. This marks a major milestone in 
the modelling of radio emission from extensive air showers. 

2. Modern modelling approaches 

In this section, we provide an (incomplete) overview of mod- 
ern modelling approaches for radio emission from extensive air 
showers and demonstrate how they can be grouped into mod- 
els predicting unipolar pulses and models predicting bipolar 
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Figure 1 : Time pulses (left) and frequency spectra (right) calculated with the large impact parameter analytical radio emission model of Gousset et al. U for the 
acceleration component (top, unipolar) and the current component (bottom, bipolar). 



pulses. For models we could not mention here, we kindly refer 
the reader to the review by Huege QJJ] . 

2.1. Models with unipolar pulses 

In 2003, Falcke & Gorham [9] motivated that the radio emis- 
sion from extensive air showers could be described as "geosyn- 
chrotron radiation" arising from the acceleration of air shower 
electrons and positrons in the earth's magnetic field. This ap- 
proach was afterwards followed by a number of authors with 
implementations in both the frequency-domain (see section 
I2.21 i and time-domain. In principle, calculations in both the 
frequency-domain and time-domain can be considered equiv- 
alent. As the radio signal, however, is very localized in time, a 
time-domain calculation has a number of technical advantages. 
Also, the treatment of retardation effects is simpler in the time- 
domain than in the frequency domain. 

The first time-domain implementation of the geosynchrotron 
approach, for the emission from the air shower maximum only, 
was made by Suprun et al. ifioll . In contrast to the frequency- 
domain calculation carried out in parallel by Huege & Falcke 
1 1 111 , the frequency spectra predicted by the Suprun et al. calcu- 
lation levelled off at very low frequencies and the pulses were 
unipolar (see Figs. 4 and 5 in lllOlO . 

In the following years, much effort was spent on developing 
a detailed time-domain implementation of the geosynchrotron 
model by Huege & Falcke and later Huege, Ulrich & Engel 
in a series of publications JI2, 13, 3. The first implementa- 
tion (REAS1, [12]) was based on a parameterized air shower 
model. Later, a transition was made to realistic air shower char- 
acteristics derived on a per- shower basis from CORSIKA [ 15] 
simulations (REAS2, 11411 '). A common feature of all of these 
time-domain calculations was the prediction of unipolar pulses 
with frequency spectra leveling off at very low frequencies, as 
seen in Fig.|2]for the case of REAS2. 



The geosynchrotron approach was also implemented in the 
AIRES air shower radio emission code II 1611 by DuVernois et 
al. iTTvh - This implementation exhibited the same characteristic 
unipolar pulses and frequency spectra leveling off at very low 
frequencies. However, the pulses predicted by this code were 
always a factor of 10-20 higher in electric field amplitude than 
those of REAS2, in spite of the near-identical implementation 
of the radiation physics. 

A recent modelling effort by Chauvin et al. Il8ll is based 
on an analytic approximation of the radio emission from air 
showers, described in the time-domain. Its results look strik- 
ingly similar to those of the various geosynchrotron implemen- 
tations. In particular, the model also predicts unipolar pulses 
(see Fig. [3} and thus frequency spectra leveling off at very low 
frequencies. Amplitude-wise, its predictions seem to agree with 
the AIRES-based implementation of the geosynchrotron model 
and are thus again a factor of 10-20 higher than those predicted 
byREAS2|fri. 



2.2. Models with bipolar pulses 

Motivated by the initial article of Falcke & Gorham, a 
frequency-domain calculation of geosynchrotron radiation was 
carried out by Huege & Falcke in As the calculation was 
based on synchrotron spectra of individual particles, the result- 
ing integrated frequency spectra showed a fall-off to zero at 
very low frequencies (cf. Figs. 6 and 9 in UHn . The frequency- 
domain implementation of the geosynchrotron model therefore 
falls in the group of models with bipolar pulses. Due to its an- 
alytic nature, however, the model had to employ a number of 
approximations and simplifications. Further modelling efforts 
therefore concentrated on an implementation of the geosyn- 
chrotron model in the time-domain Monte Carlo code REAS, 
cf. section |2~T1 
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Figure 2: Unipolar radio pulses (left) of a 10 18 eV air shower with 30° zenith angle simulated with REAS2 fl4ll and the corresponding frequency spectra (right). 




Figure 3: Unipolar pulses calculated with the ATM model fUjll . 



Another approach at modelling the radio emission from air 
showers was taken by Engel, Kalmykov & Konstantinov lfl9ll . 
They implemented a frequency-domain Monte Carlo descrip- 
tion in a variant of the EGS code [20], modelling the radiation 
from straight line-segments and the starts and ends of trajecto- 
ries in the Fraunhofer limit. The frequency spectra predicted by 
this approach showed a drop to zero at very low frequencies, 
visible in Fig. [4] Again, therefore, a frequency-domain calcula- 
tion falls into the group of models with bipolar pulses. 

With their macroscopic geomagnetic radiation model 
(MGMR, 0, HI), Scholten, Werner & Rusydi developed a 
macroscopic approach which includes the effects of the induced 
transverse current, induced dipole moments, charge excess and 
the influence of a non-unity refractive index. The dominant 
emission contribution comes from the time-varying transverse 
currents propagating through the atmosphere with nearly the 
speed of light. In contrast to the models following the geosyn- 
chrotron approach, MGMR does not treat the emission from 
individual particles, but stresses that coherent emission arises 
from the motion of charges averaged over appropriately chosen 
distances in the EAS. Consequently, in a simple approximation 
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Figure 4: Frequency spectra calculated with the EGSnrc-based radio emission 
code fl^l . The spectra drop to zero at very low frequencies, corresponding to 
bipolar pulses in the time-domain. 



where all charges are concentrated in a point 01, an analytic 
formula can be derived to relate the time-structure of the radio 
pulses directly to the longitudinal evolution of the air shower. 

Although the MGMR model describes the radio emission in 
the time-domain, it predicts bipolar pulses and frequency spec- 
tra falling to zero at very low frequencies (see Fig. [5]). It was 
this result that sparked the discussions about discrepancies be- 
tween models predicting bipolar and unipolar pulses. 

3. The missing radiation component 

As discussed in the introduction, it became clear from gen- 
eral arguments that the zero-frequency component of the ra- 
dio emission should vanish and the pulses could therefore 
not be unipolar. Only recently, however, it was identified what 
caused the discrepancies between the models predicting unipo- 
lar pulses and those predicting bipolar pulses. The reason for 
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the discrepancies is that the models with unipolar pulses miss 
a radiation contribution. All of these models start from the 
Lienard-Wiechert description of the electric field produced by 
a single moving charged particle, 



E(x, t) 



2 (1 -$-ft) 3 R 2 



+ - 



ft X [(ft - P) x /3] 



(1 -p-itfR 



(1) 



in which they only take into account the continuous acceleration 
of particles in the magnetic field. 

To calculate the radio emission from an air shower, one has to 
calculate the radio emission from an ensemble of N relativistic 
charged particles. The number N, however, varies over the air 
shower evolution and thus has a time-dependence N(t). The 
models predicting unipolar pulses effectively integrate the total 
radio emission in the air shower as: 



E tot (x, t) = N(t)E(£, t) 



(2) 



This, however, is not correct. Although this equation seems to 
explicitly take into account the fact that the number of relativis- 
tic charged particles is changing with time, a radiation com- 
ponent associated with this variation of the number of charges 
is neglected. This becomes obvious when recalling that the 
electric field equation ([TJ itself is derived from the Lienard- 
Wiechert potentials 



(D(f, f) = 
X(x, t) = 



(1 -/3-ft)R 

(1-0- ft) 



E(x, t) = -V<D(f, f) - A(x, t). 



(3) 



(4) 



The fact that not only one charged particle is radiating, but 
that the number of charged particles N(t) is changing as a func- 
tion of time has to be taken into account in the calculation of 
A(x, t) already. The time-derivative applied to calculate f) 
then leads to additional radiation terms appearing because of 
the time-dependence of N(t). 

The most straight-forward way to fix the models suffering 
from this problem is thus to not calculate the electric fields 
from single particles and integrate over them, but to calculate 
the integrated vector potential associated with the ensemble of 
charged particles A lot (x, t) and only as a last step calculate the 
electric field E tol (x, t) via equation (0]l. This approach has been 
chosen in the MGMR model |7]. Another option to take into 
account the missing radiation component has been chosen in 
REAS3 IS la] and will be discussed in the following section. 



4. Endpoint treatment in REAS3 

The radio emission component missing in REAS1/2 can be 
incorporated by considering the radiation associated with the 
"endpoints" of the individual particle tracks in the Monte Carlo 
summation. This has been realized in REAS3 In this 

implementation, the radio emission from charged particles be- 
ing deflected in the earth's magnetic field is described as ra- 
diation arising from "kinks" joined by straight line segments. 
Most importantly, however, there is additional emission at the 
starts and ends of the particle tracks associated with the "in- 
stantaneous" acceleration of particles from rest to their velocity 
near the speed of light, and the deceleration of particles from 
their velocity to rest, respectively (cf. Fig. [6}. If the number of 
charged particles at a given atmospheric depth changes, there 
will be a different number of "start-points" than "end-points", 
leading to a net contribution caused by the varying number of 
charges. This effectively takes into account the emission contri- 
bution missing in earlier implementations, and leads to bipolar 
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pulses and frequency spectra falling to zero at very low frequen- 
cies (cf. Figs|7]and[8]l. 

A very similar approach has been chosen in a recent Monte 
Carlo code for the calculation of radio emission from showers 



in dense media by Alvarez-Muniz, Romero-Wolf & Zas [ 22] . 
In general, this "endpoint" approach has been found to be a 
powerful way of describing radiation from arbitrary ensembles 
of moving particles, as has been detailed in an article by James 
et al. 12311 . The particular strength of this approach is its uni- 
versality, which does not require any differentiation between 
emission mechanisms such as synchrotron radiation or charge 
excess emission. The revised implementation in REAS3 thus 
incorporates not only synchrotron emission, but describes the 
complete radio emission arising from the underlying charged 
particle motion. 




v=0 



Figure 6: In REAS3 d, radio emission is calculated from the acceleration 
occurring in "kinks" joined by straight line segments. At the start and the end of 
a particle trajectory, additional radiation is generated due to the "instantaneous" 
acceleration of particles from rest to v » c and vice- versa. 



5. A comparison between REAS3 and MGMR 

Now that the discrepancies between the geosynchrotron ap- 
proach and the macroscopic geomagnetic radiation model have 
been resolved, we proceed to make a detailed comparison be- 
tween REAS3 and MGMR radio emission simulations. 

5.7. Methodology 

To carry out a direct comparison, we have defined a set of 
prototype air showers with a fixed geometry, a specific parti- 
cle energy and a fixed set of observer locations. Shower-to- 
shower fluctuations have been removed by choosing one typi- 
cal air shower from a set of CORSIKA simulations which was 
consecutively used as the basis for both the REAS3 and MGMR 
simulations. The magnetic field and observer altitude were con- 
figured with values appropriate to the southern site of the Pierre 
Auger Observatory. It should be noted that the REAS3 simula- 
tions shown here exhibit a slightly elevated level of numerical 
noise. This applies mostly to near-vertical air showers and will 
be improved in the release version of the code. In the MGMR 
model, numerical noise is less of an issue since the motion of 
charges is averaged in the initial stage of the calculation of co- 
herent radiation, and afterwards the field is calculated. 



5.2. Vertical 10 17 eV showers 

We first compare radio pulses emitted by a vertical 10 17 eV 
air shower in Fig. [7] Both the bipolar pulse shapes and the 
amplitudes agree remarkably well, with slight deviations close 
to the shower axis. Likewise, the agreement between the fre- 
quency spectra depicted in Fig.[8]is very good with only slight 
discrepancies at small distances from the shower axis. A look 

vertical, 10 17 eV, Argentina B-field, 1400 m a.s.l. 
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Figure 8: Comparison of the absolute value frequency spectra for vertical show- 
ers with a primary energy of 10 17 eV for REAS3 (thick lines) and MGMR (thin 
lines). 

at the lateral distributions of the unfiltered (i.e., unlimited band- 
width) radio pulse amplitudes in Fig.|9]shows that both models 
exhibit a noticeable east-west asymmetry. This is caused by 
the variation of the net charge excess in the air shower, caus- 
ing a radiation contribution with a radial polarization signature 
|H H |H 13 . It is thus clear that the radio emission from exten- 
sive air showers cannot follow a pure geomagnetic v x § polar- 
ization. The same signature can be seen in contour plots of the 
60 MHz spectral emission component shown in Fig. [TO] which 
would exhibit a strict east-west polarization for pure geomag- 
netic emission. The presence of such a deviation from purely 
geomagnetic emission is currently being investigated with ex- 
perimental data, as presented for example in 1 25] . 

5.3. No magnetic field 

If the magnetic field is switched off completely or is cho- 
sen parallel to the air shower axis in both REAS3 and MGMR, 
only the pure radio emission from the varying charge excess re- 
mains. It shows the expected radial polarization characteristics 
and is extremely similar for both models as shown in Fig. QT| 
This once more clearly illustrates that the "endpoint" approach 
implemented in REAS3 universally predicts the complete radio 
emission, even the components not associated with magnetic 
field effects, in agreement with the MGMR predictions. 

5.4. Inclined showers 

In Fig. [T2] radio pulses from a 10 17 eV air shower with 50° 
zenith angle are compared between REAS3 and MGMR. A no- 
ticeable discrepancy between the two models arises close to the 
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Figure 7: Comparison of the western polarization unlimited bandwidth radio pulses for vertical showers with a primary energy of 10 17 eV for REAS3 (left) and 
MGMR (right). 
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Figure 9: Comparison of the absolute field strength unlimited bandwidth pulse amplitude lateral distribution for vertical showers with a primary energy of 10 17 eV 
for REAS3 (left) and MGMR (right). 



shower core, where RE AS 3 predicts significantly larger am- 
plitudes than MGMR. This is not surprising, as close to the 
shower axis, details of the air shower model become increas- 
ingly important. The MGMR model currently uses a somewhat 
simplified parameterization for the thickness of the air shower 
"pancake" and neglects the lateral distribution of particles. For 
inclined air showers, this has a stronger influence than for near- 
vertical showers, as identical ground distances correspond to 
smaller effective axis distances. The deviations of the two mod- 
els near the shower core can also be seen in the unfiltered (un- 
limited bandwidth) lateral distributions of the radio pulse am- 
plitudes shown in Fig. [13] 

6. Conclusions 

After many years of research, the modelling of radio emis- 
sion from cosmic ray air showers has achieved a true break- 
through. Discrepancies between models predicting unipolar 
pulses and models predicting bipolar pulses have been recon- 



ciled. They were caused by a missing radio emission contribu- 
tion arising from the variation of the number of charged parti- 
cles during the air shower evolution. This emission component 
can be incorporated in existing models by calculating the vector 
potential and deriving the electric field only as a last step in the 
calculation, as has always been the case in the MGMR model. 
Alternatively, a powerful and universal approach using "end- 
points" can be used to take the missing radiation into account. 
The latter approach has been followed in REAS3. Compar- 
isons between REAS3 and the MGMR model confirm that the 
agreement between these two completely independent models 
based on very different approaches is indeed good. This is a ma- 
jor milestone, since previously different models predicted both 
qualitatively and quantitatively vastly different results. Remain- 
ing discrepancies in the radio emission amplitudes close to the 
shower core predicted by RE AS 3 and MGMR can be attributed 
to differences in the underlying air shower model, which is 
parameterized in MGMR while it is coming from CORSIKA- 
derived histograms in case of REAS3. 
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Figure 1 1 : Pulses for vertical 10 eV in absence of magnetic field or magnetic field parallel to the shower axis for REAS3 (left) and MGMR (right). 
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